Nanoscale Zeeman localization of charge carriers in diluted magnetic 
semiconductor-permalloy hybrids 
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We investigate the possibility of charge carrier localization in magnetic semiconductors due to 
the presence of a highly inhomogeneous external magnetic field. As an example, we study in de- 
tail the properties of a magnetic semiconductor-permalloy disk hybrid system. We find that the 
giant Zeeman response of the magnetic semiconductor in conjunction with the highly non-uniform 
magnetic field created by the vortex state of a permalloy disk can lead to Zeeman localized states 
at the interface of the two materials. These trapped states are chiral, with chirality controlled by 
the orientation of the core magnetization of the permalloy disk. We calculate the energy spectrum 
and the eigenstates of these Zeeman localized states, and discuss their experimental signatures in 
spectroscopic probes. 

PACS numbers: 73.20.-r,75.50.Pp,75.75.+a 



Diluted magnetic semiconductors (DMS) based on III- 
V alloys doped with Mn have attracted a lot of interest 
recently, due to their relatively high Curie temperatures 
T c (110 K for Gao.95Mno.o5As)|l|, below which they ex- 
hibit ferromagnetic order. In the ferromagnetic state, 
the charge carriers are spin-polarized, making these ma- 
terials ideal sources of spin-polarized currents. To date, 
most suggested applications involving DMS are based on 
this property and therefore are restricted to operation in 
a range of temperatures T <C T c . On the other hand, 
in the paramagnetic phase T > T c both the IIIi-zMn^V 
and the more established lli- x Mn x Vl DMSs (which have 
even lower critical temperatures) show giant Zeeman re- 
sponse to external magnetic fields. In our opinion, this 
can lead to interesting applications in the paramagnetic 
state, and consequently at rather elevated temperatures. 

Convincing experimental evidence for Zeeman splitting 
in the range of 30 meV for external fields of a few Tesla 
is provided by photoluminesce spectroscopy studies 0- 
Even relatively small external magnetic fields of 0.1 — 0.5 
T can easily lead to a 15 meV splitting of electronic en- 
ergy levels 3]. Comparing this to the vacuum Zeeman 
splitting of ~ 0.06 meV (for B = 0.5T) suggests that the 
effective gyromagnetic ratio of charge carriers in diluted 
magnetic semiconductors is g > 500. The origin of this 
hugely enhanced Zeeman effect is attributed [3j to the 
strong magnetic coupling J2i Js P ~d{r — Ri)s- Si between 
the spin s of the charge carrier and the spins Si of the 
Mn located at R4. In the paramagnetic state, a small 
magnetic field B induces a magnetization (Si) ~ x& 
of the Mn spins, resulting in an effective Zeeman-like 
s ■ B coupling between the charge carrier spin and the 
magnetic field, in addition to the regular Zeeman cou- 
pling —go^BS- B present in non-magnetic semiconduc- 
tors. The scale of this additional coupling is set by the 
large exchange energy J sp -d and results in a large effec- 
tive <?-value. This also implies that g(T) has a strong 



T-dependence through the magnetic susceptibility, and 
therefore can be tuned over a large range of values. This 
scenario is strongly supported by magneto-optical ab- 
sorption measurements |4| of the Zeeman splitting at the 
band edge, which clearly exhibits a Brillouin-type depen- 
dence on the magnetic field. 

The presence of giant Zeeman response in DMS implies 
that a moderate external magnetic field with a strong 
spatial variation on nanometer scale can be a very effec- 
tive confining agent for spin-polarized charge carriers in 
these systems. Provided that such highly inhomogeneous 
external fields can be created and controlled adequately 
in a DMS, Zeeman- induced localization presents a new 
route for manipulation of spin-polarized charge carriers 
at relatively high temperatures. 

Non-uniform magnetic fields with nanoscale spatial 
variations are known to appear in a variety of sys- 
tems, such as the Abrikosov flux lattice []| and arrays 
of nanoscale holes in superconducting films Q. How- 
ever, one of the most promising possibilities, which we 
investigate in this Letter, is provided by the magnetic 
vortex state Q of nanoscale magnetic disks of ferromag- 
netic permalloy NigoFe2o, Co or Fe. In the remanent 
state of such nanomagnets, the local magnetization near 
the perimeter of fairly thin disks has an in-plane vortex- 
like arrangement. This is energetically more favorable 
than a single ferromagnetic domain, since the exchange 
energy lost due to the gradual vortex rotation is more 
than compensated by the cancellation of the total dipole 
energy. This pattern is maintained for almost the en- 
tire volume of the disk. However, near the center of the 
disk, exchange interaction wins over dipole-dipole inter- 
action and shape anisotropy, and the local magnetization 
is forced out of the plane of the disk. What is remarkable 
about this topological singularity in the magnetization is 
the extremely short length scales and high fields involved: 
Very recent experimental investigations indicate that the 
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radius of the magnetic core is about 30nm in permalloy 
disks and 10 nm in Fe disks |l0| . with maximum 

field values at the core in the 0.5-1.0 Tesla range. In 
this paper we show that such highly inhomogeneous mag- 
netic fields provides a very effective localization agent for 
charge carriers in a system with large Zeeman effect. 

We investigate the properties of magnetic disks for 
which the height d is small compared to radius R. In 
this limit, the disks exhibit vortex magnetization of the 
following general type [lll | 



M(r f ) = M <f> (r)e <j ,+M z (r)e z 



(1) 



where cylindrical coordinates r — (r, <f>, z) are used. [Note 
that r 7^ \ f\]. Using Maxwell's equations [l2j, we find the 
magnetic field created by a magnetization of the type 
described in Eq. (Q. For the setup shown in Fig. ^ 
the magnetic field created in the DMS layer (z > 0) is 
B{r) = b(r) — b(r + de z ), where, in cylindrical coordinates 
f = {r,(f>,z): 



b r{r, z) = 



dppM z (p) 
o \f(r + p) 2 + z : 



{K(f(r,p,z)) 




r (nm) 



r (nm) 



FIG. 1: B z (r, z) (left) and B r (r, z) (right) in units of poMo « 
1 T, plotted as a function of r for several values of z, for a 
disk of height d = 50 nm and a core radius p c — 30 nm. The 
discontinuity in the slope of B z (r, 0) at r = p c is a consequence 
of the discontinuity of dM z {r) / dr at r = p c [see Eq. £J]. 



-E(f(r,p,z))-[p 2 + z 2 



b<f,(r, z) = 



(3) 



and 



Hoz f pc dppM z (p) E{f{r,p,z)) 

Z[r,Z '~ - Jo {r-pY+Z*- Vir + p) 2 +z 2 W 

while the corresponding magnetic vector A(r) — a(r) — 
a(f+ de z ), in the Coulomb gauge V • A = 0, is given by 



and 



a r (r, z) = a z (r, z) = 
p z f Pc dppM z {p) 



(r - p)U 



Arp 



{r + pf 



Here, f(r,p,z) = y/(4rp)/[(r + p) 2 
tic functions K(k), E(k) and II(V, k) are defined in Ref. 

The fact that only M z (r) enters these equations 
is expected, since the field lines induced by M<j,{r) are 
closed, and therefore this component does not contribute 
to the magnetic field outside the magnetic disk. Recent 
micromagnetic simulations have shown [ToL llif that the 
magnetization inside the disk is well fitted by the simple 
parameterization [ill H5I fl6| : 



[(r - p) 2 + z 2 ]^ 1 } (2) an d \M<p(r)\ = \J Mq — M 2 (r) . Such a structure is also 

known as a meron configuration in quantum Hall systems 
|l7j . In permalloy disks, the saturation magnetization is 
PqMq = 1.06 T and the core radius p c w 30 nm, if the 
disk radius R ^> p c . [The core radius p c decreases with 
decreasing d]. The core magnetization M z has been ob- 
served experimentally in permalloy disks of typical height 
d = 50 nm and radii R = 0.1 — 1 (im Q as well as Fe 
disks with d = 9 nm and R = 200 - 500 nm 

Using Eqs. (121 II) and JJJ), we find the magnetic field 
created by the vortex inside the diluted magnetic semi- 
conductor to have a typical profile as shown in Fig. QJ. 
Inside the DMS layer, B z (r, z) is largest at r = z = 0, and 
decreases with increasing r and z. B r (r = 0, z) = as 
expected for a dipole-like field. However, B r (r, z) reaches 
a maximum value which is roughly half of the maximum 
B z (r, z) value for the same z, at a distance r of the or- 
der of p c . The giant Zeeman effect present in the DMS 
systems leads to the possibility of trapping carrier states 
in the strong magnetic field near the disk surface. 

The eigenstates of a charge carrier placed in such a 
magnetic field are given by the Schrodinger equation: 



(5) 



(r + p)y/{r + p) 2 + z 2 



f(r,p,z) \+(r + p)K(f(r,p,z)) 



(6) 

2 ] and the ellip- 
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2m 



p + qA{r) 



Em (8) 



\M z (r)\/M = 







r 2 )/(p 2 + r 2 ) ,r<p c 
,r>p c 



(7) 



where, as already mentioned, the large ^-factor effec- 
tively accounts for the magnetic exchange of charge carri- 
ers with locally magnetized Mn spins. The Hamiltonian 
given in Eq. (JHJ neglects electron-electron interactions. To 
first order, this is justified if the carrier concentration is 
either large enough that screening is very effective, or 
for low electron concentrations where only a small num- 
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ber of electrons might be trapped below each disk. In 
the latter case, Eq. (JSJ) should accurately describe the 
strongest-bound state; for higher energy bound states, 
one should also include the screening provided by elec- 
trons occupying inner shells. We omit this complication 
here. Eq. @ also assumes that the Mn spins create a 
smooth polarizing field for the electron, ignoring possi- 
ble local fluctuations due to the fact that Mn spins are 
distributed randomly in the DMS. This so-called virtual 
crystal approximation has already been used extensively 
to describe DMS systems. 

Given the symmetries of the magnetic and vector fields 
[Eqs. 1)216(1 ]. Eq. (JHJ) has solutions of the following general 
form: 



4> m {r, <j>, z) = exp(im<, 



^\r,z) 
Ik 71 ' (r, z) exp [i<j>\ 



(9) 



where the angular momentum m is an integer. The ex- 
tra phase exp [iip] in the spin-down component has two 
significant consequences: (1) a z is not a good quan- 
tum number. [Besides the (s 2 ) ^ part, the expec- 
tation value of the spin acquires a radial contribution 
(s x +is y ) ~ exp (i4>), for all m. This is clearly due to the 
existence of a radial part B r (r 1 z) in the magnetic field], 
and (2) the (usually expected) degeneracy between states 
with ±m is now lifted. The reason is that the presence 
of the magnetic field breaks the time-reversal symmetry 
responsible for this degeneracy. This has interesting con- 
sequences, as described below. 

The general form of the wave function given above 
allows us to draw immediate conclusions regarding se- 
lection rules governing the electromagnetic response of 
the electronic states trapped in the DMS by the vor- 
tex magnetic field. Assume that the system is ex- 
posed to monochromatic radiation, and let Ai,(r) and 
= V x Al{t) be the vector potential, respectively 
magnetic field associated with it. Then the Schrodinger 
equation JHJ acquires three extra terms, proportional to 
a ■ Bl(t), (p + qA) ■ Ai(r) and to A 2 L (r), in order of de- 
creasing magnitude. Assume that the beam propagates 
along the z-axis and is circularly polarized, B^ (z) ~ 
[e x ± iey) exp (ikz). Using Eq. it is straightforward 
to find the selection rules (m'\a ■ |m) ~ £m',m±i- 
(The term (p + qA) ■ A^if) obeys the same selection rules, 
while the A 2 L (r) term induces two-photon processes but 
with a vanishingly small probability 18]). This selection 
rule implies that the absorption of a right (left) circularly 
polarized photon excites the electron to a level with an m 
increased (decreased) by one unit. If levels with ±m are 
no longer degenerate, the system will interact differently 
with photons of different circular polarizations. 

Let us now use a simplified model to demonstrate the 
lifting of the ±m degeneracy. First, the terms involving 
the vector field A(r) are removed from Eq. JSJ. This is 
justified since they are vanishingly small compared to the 



Zeeman term, due to the supplementary enhancement of 
the later by the interaction with the Mn spins. Second, 
since we are interested in the most strongly-bound states, 
which are likely to be localized at small (r, z) values, we 
use Taylor series for the magnetic fields in this region, 
given (up to quadratic terms) by: 



B r (r,z) 
B z (r, z) 
where the functions 
A{z) = b x - b 2 z + 



MoM 5 ( — ) — 

\PcJ Pc 



PoM 



A | — 

Pc 



6a 



26, z 2 



and B(z) 



2 



2hz 



have been introduced for later convenience. These fields 
continue to satisfy the condition V • B = 0, as neces- 
sary. The coefficients b\ , 62 and 63 are complicated func- 
tions of d/p c . In the limit d/p c — > 00, we find b\ = 1/2, 
62 = (7r + 2)/4 and 63 = 1. These asymptotic values turn 
out to give reasonable approximations for the typical ex- 
perimental parameters d — 50 nm and p c — 10 — 30 nm, 
with relative errors less than 10% for d/p c = 5/3 and 
decreasing very fast to zero for larger ratios [19]. As a 
result, the asymptotic values will be used in the remain- 
der of this paper. 

With these simplifications, we look for eigenfunctions 
[see Eq. (0] of the following form: 



4>^\r,z) 



= ax (zip, 

0,2 (z/Pc) 



m 



exp 



|m+l| 



exp 



b 2 p 2 c 

v 2 
b 2 p 2 c 



Both components are regular for r — > 0. The functions 
ax(z/p c ) and ci2(z/p c ) are determined, in dimensionless 
units zj p c — * z, by the set of coupled equations: 



4(|m| + l) 



b 2 



b 4 



dz 2 



-aA(z) 



s-1 

a\ H — aB(z)a2 = eai 

(10) 



a\(z) 



-aB(z)a,2(z) = 



(11) 



A(\m\+s+l) d 2 

b 2 ~d^ +aA{z) 



3 + 1 uf \ 
a,2 aB{z)a,i — ea 2 



(12) 



( 4 \ s-1 
- {-^ + a\a 2 (z) + ^-aB(z)a 1 (z) = (13) 

Here, s — 1 if m > and s = — 1 if m < 0. We 
define the energy unit Eq — h 2 /(2mp^.) and use e = 
E/Eq, while 2a = gp,BPoM /E is the ratio between 
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the maximum Zeeman energy and this energy unit. This 
number is rather large. If we use p c — 30 nm and /io-Mo = 
1.06 T (typical values for permalloy disks), to — 0.5m 
(band value for heavy hole mass in GaAs) and use the 
estimated effective gyromagnetic factor to be g = 500, 
we find |a| ~ 175. 

If s = +1, we see from Eq. (|13|l that there is a non- 
trivial solution if and only if 4/6 4 = —a, which is possible 
if a < 0. This is the case when Mq < 0, i.e. the z-axis 
disk magnetization M z points away from the DMS sur- 
face [13. With b 2 = vWH, Eq. (TTJ has the solution 
ai(z) = 7}B(z)a,2(z). Using this in Eq. I|12[l we find 
a 2 (z) = cxp (— kz), where k > is linked to the eigenen- 
ergy through e m = 2y/\tT\(m + 2) + a (6 ~ 7r ^ 1 8 ° +7r) - k 2 . 
Finally, Eq. (fTUfl must be satisfied up to powers of z° 
(the magnetic fields, and therefore the functions a\ and 
aa, are accurate only up to z 2 , meaning that the second 
derivative of a\ is accurate only up to z°). This gives 
k = -Q.321- N /Ia[+0.127|a| (with k > satisfied for any 
\a\ > 1.56). In other words, up to a large negative con- 
stant the spectrum e m = 2^/|a|m is similar to that of a 
harmonic oscillator (s = +1 — > m > 0) |2l|. For typical 
values for the vortex core radius and effective mass (con- 
stants which set the energy scale E as given above), the 
level spacing is in the range of 3 — 6 meV, which should 
be accessible by most spectroscopic tools. 

To summarize, for M z < we find only solutions with 
s = +1 or m > 0, proving that the ±m degeneracy is 
indeed lifted in this simplified case. For M z > 0, the 
reverse is true; we find only solutions with m < 0. How- 
ever, these simple solutions only hold for small values of 
to, where the wave-functions are localized at small r and 
z and the Taylor series for the magnetic fields are valid. 
To find the true spectrum, one must integrate Eq. 
numerically for the full expression of the magnetic fields; 
this work is in progress and the results will be reported 
elsewhere We expect that the degeneracy between 
±m eigenstates is lifted in the general case as well. As al- 
ready discussed, this means that the system will interact 
differently with photons of different circular polarization, 
depending also on the orientation (sign of) the magnetic 
disk core magnetization M z . Such a property may be 
useful in designing spintronic devices : the memory bit 
(M z up or down) controls the properties of the trapped 
electronic states, and therefore the behavior of the whole 
device. Other geometries, such as lines of ordered disks, 
could be used to trap and transport electric currents with 
similar non-trivial spin properties. 

In conclusion, we suggest a different route to creat- 
ing spintronic devices that can operate at room tem- 
perature, by combining the giant Zeeman effect of di- 
luted magnetic semiconductors in the paramagnetic state 
with the highly inhomogcncous magnetic field created by 
nanoscale permalloy disks, or other nano-patterned mag- 
netic material structures. In particular, we have shown 
that electronic states can be trapped near the surface of 



a magnetic disk. Their properties can be suitably tai- 
lored using various diluted magnetic semiconductors and 
various types of magnetic disks, as well as by varying 
the temperature (since the effective g-factor has signifi- 
cant temperature dependence). Since time-reversal sym- 
metry is broken in the presence of the magnetic field, 
the system will interact non-trivially with any other sys- 
tem which has a definite chirality, such as circularly po- 
larized light. To our knowledge, this is the first time 
that Zecman-induced localization has been suggested and 
demonstrated theoretically, or that electronic states with 
such unusual spin-polarization have been derived. Ex- 
periments to test these i are currently in progress |22^ . 
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